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The transition from linear to nonlinear dynamical elasticity in rocks is of considerable interest in 
seismic wave propagation as well as in understanding the basic dynamical processes in consolidated 
granular materials. We have carried out a careful experimental investigation of this transition 
for Berea and Fontainebleau sandstones. Below a well-characterized strain, the materials behave 
linearly, transitioning beyond that point to a nonlinear behavior which can be accurately captured by 
a simple macroscopic dynamical model. At even higher strains, effects due to a driven nonequilibrium 
state, and relaxation from it, complicate the characterization of the nonlinear behavior. 
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Rocks possess a variety of remarkable nonlinear elas- 
tic properties including hysterisis with end-point mem- 
ory PljVariation of attenuation and sound velocity with 
strain Q , strong dependence of elastic and loss constants 
on pressure, humidity, and pore fluids yl, long-time re- 
laxation phenomena ('slow dynamics') and nontriv- 
ial variation of resonance frequency with strain Q. 
Significantly, materials as diverse as sintered ceramics 
and damaged steels are now known to display similar 
effects 0. Thus rocks may be viewed as representative 
members of a class of fascinating, but poorly understood, 
nonlinear elastic materials: Fundamental questions still 
to be resolved relate not only to the underlying causes 
of the nonlinear phenomena but also to the conditions 
under which they occur. 

In this Letter, we focus on delineating two strain 
thresholds, one below which the rocks behave effec- 
tively as linear elastic materials, e^, the other beyond 
which memory and conditioning effects occur, ejvf, and 
the dynamic elastic behavior straddling the region of 
these thresholds. While in Ref. 2\ it was argued that 
El ~ 10~^ (albeit with some uncertainty), more recent 
data 0, IE H have been used to support an extension of 
the nonlinear region to substantially lower strains; doubt 
has been cast even on the very existence of a thresh- 
old 1^. In addition, results from resonant bar experi- 
ments [E have been interpreted to exhibit a 'nonclas- 
sical' frequency and loss dependence on the drive ampli- 
tude, i.e., frequency and Q softening linearly with drive 
amplitude rather than quadratically as predicted by Lan- 
dau theory 0, even at strains as small as 10~^. (The 
importance of ej\/ in interpreting resonant bar data is 
emphasized below.) 

We have carried out a new set of well-characterized ex- 
periments, over a wide dynamic range, to unambiguously 
settle these questions: While longitudinal resonant exci- 
tation of bars is a classic measurement technique |lOj . 
rock samples require substantial care in terms of con- 
trolling the temperature and humidity and characteriz- 



ing possible systematic effects, especially those due to 
conditioning of the sample by the external drive |^ . 

Our major conclusions are as follows. For Berea 
and Fontainebleau sandstone samples, below a threshold 
strain cl ^ 10~^ — 10~^ (lower end for Fontainebleau, up- 
per end for Berea), there was no discernible dependence 
of the resonance frequency on the strain - the materi- 
als behaved linearly to better than 1 part in 10''. For 
£>£!,, the materials displayed a reversible softening of 
the resonance frequency with strain, in excellent qualita- 
tive agreement with the quadratic prediction of classical 
nonlinear theory |^ up to a point where memory and 
conditioning effects became apparent (e ea/)- In de- 
tailing and characterizing the onset and nature of the 
nonlinearity, our results very substantially improve on 
previous work We show below that, up to the condi- 
tioning threshold, the dynamical behavior is accurately 
captured by a phenomenological macroscopic model in- 
corporating a (softening) Duffing nonlinearity and linear 
losses. Beyond the conditioning threshold, the simultane- 
ous presence of nonlinearity and nonequilibrium dynam- 
ics complicates the characterization of dynamical behav- 
ior; in the absence of a separation of these effects, the 
data cannot be interpreted to support the existence of 
nonclassical behavior. 

Our computer-controlled resonant bar experiment uses 
cylindrical sandstone samples, 2.5 cm in diameter and 
35 cm in length. The cylinders are driven sinusoidally 
at one end by a PZT transducer with a brass backload. 
The acceleration is measured by a B&K accelerometer at 
the other end of the bar and converted to an effective 
average strain using the (known) driving frequency, /, 
via e = il/ {AirLf'^), where u is the displacement and L is 
the bar length. The finite accelerometer signal to noise 
restricts the useful strain sensitivity to ~ 10^ — 10^^ 
while the upper end of strain is limited primarily by the 
physical integrity of the experimental setup to ~ 10~^. 
Scans of the resonance peak are conducted at constant 
drive amplitude over an up/down frequency sweep with 
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FIG. 1: The average strain amplitude e as a function of drive 
frequency for Fontainebleau sandstone. The reference center 
frequency is 1155.98 Hz. The open circles are the experimen- 
tal data; the filled circles mark the peak positions. The solid 
lines are theoretical predictions from Eqn. (|5J. 



FIG. 2: Representative data for the resonant frequency shift 
Af as a function of the effective strain e for three samples. 
Acrylic is a linear material used as a control in the experi- 
ments. In this figure, the error bars are smaller than the sizes 
of the symbols. 



spot measurements at frequencies fi bracketing the bar's 
resonance frequency (Fig. While in previous exper- 
iments the temperature was actively controlled, the 
present experiments were carried out using passive ther- 
mal isolation to avoid even low-level thermal shocks. 

The thermal history of a rock sample is known - due 
to slow build-up and relaxation of internal stresses - to 
influence the effective elastic modulus, and hence the res- 
onance frequency of the bar. The timescales associated 
with this behavior can be quite long, of the order of sev- 
eral hours, hence long-term temperature stability is a 
basic necessity for resonant bar experiments, especially 
at small strain levels. With the present isolation system, 
long-term frequency stability of the samples has been ver- 
ified at ~ 0.1 Hz (corresponding to a long-term thermal 
stability of ~ 10 mK), which is close to how well the peak 
of the frequency response curve can be determined at the 
lowest levels of strain reported in this Letter. 

The basic quantity measured in these experiments is 
the resonance frequency, /o — wq/Stt, of the bar as a func- 
tion of the strain, e, defined by the peak of the resonance 
curve as measured above. The statistical analysis used 
a nonparametric Gaussian process to model data trends. 
Bayesian estimation and characterization of uncertainty 
was carried out using Markov chain Monte Carlo ITlj . 

Traditionally, the loss, as represented by the Q of the 
bar has also been reported. However, this quantity is not 
easy to define or measure precisely [the present bars have 
Q ~ 66 (Berea) and Q ~ 132 (Fontainebleau)] even at 
high values of Q We find experimentally that the 

full- width of the response curves, F (measured at ao/\/2, 
where oq is the peak amplitude; Q = ujq/T) at all drive 
levels below the conditioning threshold, is essentially in- 
variant (this behavior is also predicted by the macro- 



scopic model described below). At higher drives, as con- 
ditioning and memory effects appear, the response curves 
are not symmetric around the peak and the width be- 
comes strongly dependent on the drive amplitude. This 
behavior will be detailed elsewhere [T^ . 

The determination of /o as a function of strain is com- 
plicated by the fact that external driving can condition 
the sample and also lead to shifts in the frequency via 
slow dynamics 0|. In order to eliminate this possibility 
we implemented a 'zig-zag' strategy of systematically in- 
creasing the drive level through the up/down frequency 
sweeps, and then dropping back to the original drive am- 
plitude to verify that /o at the lowest strain value had not 
changed. Application of this method shows that there are 
no conditioning effects for e < eM — 5 x 10~^ for Berea 
and e < em — 2 x 10^^ for Fontainebleau. Up to these 
strain levels, any variation in /o as a function of drive 
amplitude is taken to define the intrinsic nonlinearity of 
the sample. 

The measured variation of /o over a relatively wide 
range is shown in Fig. [51 (A/ is the difference between the 
measured value and the starting, lowest drive level, value 
of /o). Results from a known linear material. Acrylic, are 
shown in order to provide a control reference. Within the 
error bars of the experiment, the materials are effectively 
elastically linear (no softening of /o with drive) till e ^ 
10~® — 10^^. Beyond this point, /o softens quadratically 
with increasing drive amplitude until e ~ cm (Fig. 
We have compared our new results with archival data 
from previous experiments on Berea reported in Refs. 0; 
allowing for differences in the starting values of /o, we 
find that the two datasets are in excellent agreement. 
This is quite remarkable, and encouraging from the point 
of view of sample-independence, since the samples were 
subjected to very different environmental conditions in 
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the two cases. The erroneous conclusion of nonclassical 
nonUnearity (linear dependence of /q on e) reached in 
this earlier work was due to the limited dynamic range 
in strain of the analyzed data (< 1 decade) versus the 
present measurements (3 decades): The present data are 
consistent with a zero value for the cofRcient of a linear 
fit to the softening with strain; the corresponding values 
of Refs. are an order of magnitude too large to agree 
with the measurements. 

We now describe a simple phenomenological dynami- 
cal model for driven rods that provides an excellent de- 
scription of the measured data. The model is not di- 
rectly extracted from a one-dimensional (nonlinear) wave 
equation Rather, the procedure follows statistical 

mechanics-based modeling of degrees of freedom coupled 
to dissipative channels 15,16]. The model adds a quartic 
(Duffing) softening nonlincarity to a harmonic potential 
(a more complex model for the higher strain regime is 
described in Rcf. 16] ). The equation of motion for the 
displacement is taken to be: 



, + 2^u + ^fu + = F sin(Ljt), 



(1) 



where 7 < 0. Since the displacement, u, is small, multi- 
scale perturbation theory can be used to solve Eqn. (Q) 
very accurately For the case relevant to the exper- 
iment (a; ^ S7), the solution is m = acos(a;t + 0). The 
phase (j) is of no interest here, while the relation between 
the amplitude of the oscillation, a, and the drive ampli- 
tude, F, is given by 



crfJ a^7 
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(2) 



where oj = Vl + a. It is straightforward to show that the 
peak of the response curve has the value, = F/ (2fin), 
and occurs at the drive angular frequency u>o = ^ + cq, 
with 



3f^7 



(3) 



Thus, the model predicts a quadratic softening of the 
frequency with the drive amplitude F. In addition, by 
solving Eqn. (|5J for a, and then substituting a — ao/V2, 
it is easy to show that the width of any response curve is 
an invariant, F = 2/i. 

With these results in hand, it is straightforward to de- 
termine model parameters. As the model predicts, we 
have verified that T as measured from the experimental 
curves is constant within 1% up to the strain cm, this 
immediately determining the damping coefficient fj.. The 
(linear) resonant frequency Q and the nonlincarity pa- 
rameter 7 now follow by fitting the experimental data 
for A/ as a function of the drive using Eqn. ^ (Fig.|2I). 
As apparent from Fig. ^ with these parameters fixed 
as just described, the model predictions are in excellent 
agreement with the experimental response curves. 

It has been previously claimed that the absence of fre- 
quency softening is not sufficient to rule out nonlincarity 
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FIG. 3: The resonant frequency shift A/ as a function of 
the effective strain e for Fontainebleau and Berea samples for 
€ < ea/. The solid lines represent predictions of the theo- 
retical model [Eqn. @]. Parameters for Fontainebleau are 
Q. = 7262.8 rad/s, fi = 27.5 s-\ 7 : 



and for Berea, f2 = 17375.7 rad/s, /i = 131.6 s 
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in rocks as harmonic generation may exist even in the 
absence of a discernible frequency shift l^: At least in 
the materials studied here, this does not occur. Our the- 
oretical model assumes that the fundamental mode dom- 
inates the response of the bar to external driving and no 
higher harmonics are excited via mode-coupling. With 
the parameters fixed as above, harmonic generation via 
the intrinsic nonlincarity of the model is very weak, with 
all even harmonics suppressed, and with odd harmonics 
typically 80 dB below the fundamental - lower than the 
noise floor of the experiment, and, consistent with this 
prediction, we did not observe any harmonic content in 
the signal even at the highest drive levels. 

We now turn to the behavior of /o as e is increased 
beyond e^v/, the point where conditioning begins to play 
a role. Here the frequency softening is such that /o no 
longer returns to the starting point when the drive is 
removed, but to a lower value. It is clear that in this 
regime the measured softening at a given strain cannot 
be interpreted as wholly due to an intrinsic nonlincarity; 
it will as well depend on the sweep-rate Analyzing 
a /o vs. e curve without taking this effect into account 
would introduce an unknown systematic effect exagger- 
ating the actual nonlincarity. Fortunately, experimental 
protocols can be implemented to alleviate this difficulty 
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FIG. 4: Avoiding relaxation effects: The left curve is for Berea 
maintained in a nonequilibrium steady state by continuous 
driving, while the right curve allows for a long relaxation to 
take place between the individual data points. 

as described in detail in Ref. To illustrate these is- 
sues, Fig. ^ shows response curves for Berea sandstone 
in two different regimes. The left curve represents an 
experiment where, at each frequency, the sample is al- 
lowed to come to a new steady state (10 minutes at each 
frequency), while the right curve allows for 10 minutes 
between rapidly taken individual measurements to allow 
a return to thermal equilibrium between measurements. 
In both cases, the measurement protocols are designed 
to reduce relaxation effects as much as possible, in one 



case by staying at all times in a nonequilibrium steady 
state, and in the other, by allowing a close-to-equilibrium 
return between measurements. The key point is that at 
strain levels e > em , the rock transitions to a nonequilib- 
rium state characterized by a different set of macroscopic 
parameters, as evidenced by the 5 Hz shift between the 
left and right curves. If measurements could be made 
rapidly enough so that each point on the right curve 
corresponds to an equilibrium state, then it is precisely 
these response curves that are equivalent to data taken 
at e < sm- In practice, however, this procedure is very 
difficult to carry out as it requires stringent long-term en- 
vironmental stability. Thus, all measurements to date in 
this higher strain regime are dangerously contaminated 
by conditioning; the separate effects of nonlinearity and 
relaxation in these experiments cannot be disentangled. 

To summarize, our experiments have established the 
existence of a reversible dynamic quadratic nonlinear- 
ity in the Berea and Fontainebleau sandstones up to 
a material-dependent strain threshold, cm- Below this 
threshold we find no evidence for nonclassical behavior 
as reported previously 0] ■ Frequency shifts in dynam- 
ical experiments with e > cm do not have a simple inter- 
pretation due to the existence of a driven nonequilibrium 
state with differing macroscopic parameters; because of 
the competition of material nonlinearity and condition- 
ing and relaxation effects, present experimental data can- 
not distinguish classical from nonclassical effects in this 
regime. Experiments to do so are in progress. 

We are indebted to Tim Darling, Robert Guyer, and 
Tom Shankland for stimulating discussions. 
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